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Conclusions 
— MMH property tables suitable for 
use in GFSSP were developed. 

— MMH properties were used to 
analyze an MMH thruster system. 
= Steady state results closely 
matched expected values, 


Introduction _ Algorithm 
Fluid Transient analysis is important for the design of spacecraft |— Obtain data or equations for defining fluid properties as a function of pressure and 
propulsion system to ensure structural stability of the system in the|| temperature; density, dynamic viscosity, enthalpy, entropy, specific heat. 
event of sudden closing or opening of the valve. Generalized Fluid)|_ For specialized liquids such as MMH, the required property data may not be available; 
System Simulation Program (GFSSP), a general purpose flow network) 4 5. density listed only as a function of T (incompressibility assumption). However, 


code Geveloped al NASA/MSFC Is capable of simulating pressure fluid transients such as water hammer occur due to fluid compressibility and this 
a due ‘to sudden opening or en of valve when ‘nformation.is vital Gi GEESP simulations. 


Nor mic Oro erties of real fluid are hie" eC Tor the eiiaire: , >. ne 
i oi mgt ons A , acataane — |f the speed of sound in the fluid is known at a Known temperature, then the pressure | indicating that the property 
"iA in order 7& predict dependence on density can be estimated; e.g. a value of 5655 ft/s was used for the tables are working in GFSSP. 


© = annie oc" urine 4 fluid transient. -Unfortune speed of sound in MMH [1]. - 
/- aman 4MiCc = a transi ar 4S REFPROP. GASP or |— MMH property data for density, p = p(p) was then generated for a fixed temperature, 
-|GASPAK do not provide the ri aaa Pr oneies of a valid assumption since the test was under isothermal conditions. 
Monomethylhydrazine (MMH). This work illustrates the process used |— Using equations [1] for Cp = Cp(T), u = u(T), and h = h(T) property tables were then r 
for building a customized table of properties of state iia bles | from; generated for specific heat, viscosity and enthalpy. | 
= 


= Transient results snow good 
agreement with 
experimental results with 
regard to the shape of the 
pressure history and the 
frequency response. 

" Oscillations decayed with 

time as expected. 

Amplitude shows some error 


available properties and speed of sound that is required by GFSSP for — GFFSP requires property values for ratio of specific heats, thermal conductivity and 
simulation. Good agreement was found between the simulations and|| entropy even though they are not used in the transient calculations. This data could 


measured data. This method can be adopted for modeling flow| pot be obtained from surveyed literature. Arbitrary values were generated as inputs 
networks and systems with other fluids whose properties are not| 


known in detail in order to obtain general technical insight. 


| tothe code. 


Preliminary Results Future Work 
Background | — Asimplified flow network of the test configuration was developed as shown below. __||—~ GFSSP modeling of the damping 
= Sy coetttnoiten iid a wes Sobek .,— When proper boundary conditions were applied, the steady state simulation of the effects of a helium bubble that 
— MOC (Method of Characteristics) uses speed of sound Poe giv MMAH test system in GFSSP was able to produce mass flow rates similar to those can develop in the MMH 
governing equation to model fluid transient | measured during the Leet : | | Thiet ecm, 
— Unlike MOC, GFSSP uses compressibility of fluid to model fluid — Atransient simulation of the system with an nearly instantaneous (1 ms) valve closure Further veritication of MIME 
- transients. was then executed: simulations in GFSSP with 
— = ~__, a i = With simplified model and proper GFFSP settings, the model runtime was experimental data. 
Objective: Create GFSSP Compatible property tables for : very quick compared to other common analysis techniques. \- More detailed property tables 


MMH and then use GFSSP to model fluid transients in 4 =" Results demonstrated similar frequency response to that seen in data from 
small MMH thruster that is under development. | an experimental test of the system (shown in plots below). 


and equations of state for MMH 
need to be developed to improve 
the accuracy of MMH transient 
analysis. 

= Lack of data over a range of 
pressures currently is a 
significant source of error. 


» Water hammer oscillations decayed at approximately the same rate as the 
test data. 

=» Amplitude of pressure oscillations was slightly different than expected. This 

was input using a customized routine implemented in Matlab. can be attributed to the geometry simplification used in the simulation and 


— Analysis based on Airforce Propellants Handbook data for MMH. due to the property (p) estimation with change in pressure. 
Equations of state were created using curve Tits to the data. Wee ae Sepesee eeeenre " Pressure vs Time At Valve 

— Created Matlab script to take the equations of state and create ae 2 —— 
GFSSP compatible property tables. 

— Code could be ecoptega for Ive other fluid. 


Pressure vs. Time at Valve (Test Data) 
—< Test005_1$1-0303-0_CUSTOMER.MT PT308 PSIA Fu aan 


Methodology 
— GFFSP requires accurate pressure-density relationship for 
simulations. Since data for MMH is not readily available, this data 
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(Left) Dynamic pressure data at the valve, from test , showing water hammer transients. 
[Above Left] Simplified GFFSP model of test configuration showing boundary conditions. 
[Above Right] Computed fluid transient at the valve due to near instantaneous valve closure, 1 ms. 


